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Abstract 

Using an exactly solvable pairing model Hamiltonian in the static path approxi- 
mation together with small-amplitude quantal fluctuation corrections in random 
phase approximation (SPA+RPA), we have analyzed the behaviour of canonical 
(number projected) and grandcanonical treatments of nuclear level densities as 
a function of temperature and number of particles. For small particle numbers 
at a low temperature, we find that though the grandcanonical partition function 
in SPA+RPA approach is quite close to its exact value, the small errors in its 
estimation causes significant suppression of level density obtained using number 
projected partition function. The results are also compared with the smoothed out 
exact values of level density. Within this model study, it appears that due to sad- 
dle point approximation to multiple Laplace-back transform, the grandcanonical 
treatment of level density at low temperature may be reliable only for relatively 
large number of particles. 
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The familiar Bethe formula[|IJ for the level density, because of its simplicity, 
has been widely used to perform the statistical analysis of nuclear reactions. This 
level density formula takes a simple form due to (i) the connection of grandcanon- 
ical partition function with microcanonical partition function (level density) by a 
saddle point approximation in the evaluation of the traces over the states[|3[] and 
(ii) the grandcanonical partition function itself is approximated using independent 
single particle spectrum which is further assumed to be equidistant. Recently, a 
more realistic value of the level density ||]-|| in saddle point approximation is ob- 
tained using a grandcanonical partition function in static path approximation^, ||] 
which accounts for the large-amplitude thermal fluctuations. The method of sad- 
dle point approximation can be generalized in order to treat the level densities 
that are characterized by a set of quantum numbers, provided, these quantum 
numbers are composed additively of contributions from the single particle states 
(e.g. excitation energy (E*), number of protons (N p ) or neutrons (N n ), angular 
momentum, parity etc.). For instance, consider the system consisting of one kind 
of particles(protons or neutrons). The exact level density, 

p(e, n) = y, 6 ( E - E M N - N *) CO 

A 

can also be obtained by double Laplace-back transform of the grandcanonical 
partition function Z(@,a) as 

p{E,N) = (— f j Z((3, a y E - aN d(3da (2) 

In saddle point approximation, the level density p can be obtained directly in terms 
of canonical and grandcanonical partition functions as 
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respectivley. The canonical partition function Z N {(3) is given as 

1 r+ioc 

Zn{(3) = — Z{(3,a)e- aN da (5) 
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In above equations, (3 corresponds to inverse temperature (1/T) and a = ^(3 with 
H being the chemical potential. Values of (3 and a in eqs. (||) and (^) are chosen 
such that the saddle point conditions are satisfied, i.e., 

dlnZ ^ dlnZ 

From eqs. (^|) and (0), it is clear that in general the canonical and grandcanonical 
treatments of level density would differ by the saddle point approximation to a sin- 
gle and a multiple Laplace-back transform of corresponding partition functions. A 
typical example, the grandcanonical treatment of level density for a given nucleus 
(i.e. fixed N p and N n ) as a function of excitation energy and angular momentum 
would require saddle point approximation to four Laplace-back transform |6|, H- 



However, it has been pointed out in ref.|I(| that the Bethe's formula for level 
density due to the saddle point approximation to several Laplace-back transform 
leads to a severe problem — violation of "microcanonical analyticity constraint" 
leading to a measurable differences between the microcanonical and grandcanon- 
ical treatments. Also, it can be seen in ref.[[5], that at very high excitation 
energy the assumption of equidistant single particle spectrum in Bethe's formula 
gives rise to exponential deviations relative to more realistic values. 

In this letter, we analyze the behaviour of the canonical and grandcanonical 
treatments of the level densities. For this purpose, we have used a recently devel- 



oped method of exact functional representation JT^, Q to construct grandcanoni- 
cal partition function for the interacting systems. It is clear from refs.[|T3|, |15|] that 



the exact functional representation can be made easily numerically tractable using 
SPA+RPA (static path approximation+random phase approximation) approach. 



Moreover, it is shown in refs. []T5|, [15[ [IT] that when SPA+RPA approach ap- 
plied to the grandcanonical partition function for the exactly solvable pairing and 
Lipkin model Hamiltonians, the results for E* , p and heat capacity found to be 
very close to ones obtained using exact grandcanonical partition function. This 
approach has been applied to a realistic Hamiltonian also[|HJ. In what follows, 
we describe briefly, the SPA+RPA approach to canonical (number projected) and 
grandcanonical partition functions for exactly solvable pairing model. 
We consider here the pairing Hamiltonian defined as 
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where G is the interaction strength, P = J2k>o The index k — 1, 2, ...Q with 
k = —k labels 2Vt degenerate single particle states. Equation (0) can be rewritten 
in terms of quasi-spin operators as 

H = -G(J 2 - J 2 Z + J z ) (8) 

where, J x = \{P + Pt), J y = \{P - pt) and J z = \(N - SI). Since, [P, P 
0, the eigenvalues can be given as 



H, J z 



where, 



E(J, M) = -G(J{J + 1) — M 2 + M) 



M = -(N-tl) 
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with iV being the number of particles. The degeneracy factor for a given J state 
is 
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With these ingredients (eqs. P-[Tl|), the exact level density can be obtained using 
eq. (1) as 

p(E,N)=J2g.AE-E(J)). (12) 

Here we have suppressed the M dependence of eigenvalues (see eq.(|)) since it is 
fixed for a given number of particles. We define now a smooth version as 

p(E,N)=Y / 9.A(E-E(J)) (13) 

j 

which describes the exact average level density. We adopt a simple form for the 
spreading function A(x) as given below, 

11 9 1 

A n {x) = -^e~ nx (14) 
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where the parameter n controls the width of spreading function. The above choice 

1 /2 2 

of A-function essentially represents the dominant part of A = L 3 ' e~ x as used 
in ref. [[TJ. 

It follows immediately from eqs. (|3|)-(|5|) that the level density in saddle point 
approximation can easily be obtained once the partition function is known. The 
exact canonical or grandcanonical partition function can be constructed as 

Z N {(3) = Y J 9. J e- pEj (15) 
j 

and 

Z(P,a) = J2e aN J29.je-^ (16) 

N J 

where J =\ M \, \ M \ +1, ...,0/2 or (Q - l)/2 depending on whether the N is 
even or odd. It is worth mentioning here that eq. ( |i~5p can precisely be obtained 
by substituting eq. (|T5| ) in eq. (|5|) and using "Wick rotation" (i.e., a ia). 
In SPA+RPA approach, Z(/3,a) takes the following form, 

Z(p, a ) = W r°° AdAe -#A3 +a n + ns + 2n M i+e- s W A ) (17) 
G Jo 
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where, the factor C(A) as given below 
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is due to RPA, which accounts for small-amplitude quantal fluctuation corrections 
and Q = a + § G, S = y/Q 2 + /3 2 A 2 . With C = 1, eq. flT7|) reduces to the 
partition function in SPA approach. The superscript tilde on Z in above equation 
is used to distinguish it from exact one. 

The canonical (number projected) partition function can be obtained as follows 



in a ensures that the partition function Z N {(3) will be real even if the approximate 
form for grand partition function is used. 

We analyze now the saddle point approximation to the level density p, eval- 
uated using (a) SPA+RPA to canonical (number projected), (b) SPA+RPA to 
grandcanonical, (c) exact canonical and (d) exact grandcanonical partition func- 
tions. For the sake of comparison of these results with the smoothed out exact 
value of level density (eq. [13]), we have plotted in figures [1] and |] the ratio p/p as 
a function of temperature for N = 5 and 10, respectively. To facilitate the further 
discussions, we denote p for the cases (a), (b), (c) and (d) by p c Tpa , p 9 rpa , p c ex and 
p^ x , respectively. 

It can be seen from figures [1] and ^] that for all the cases, p/p is quite different 
from unity at low temperature (T « 0.25 MeV) but with increase in temperature 
it approaches unity. The negative slope in p/p vs T curves indicates that the 
saddle point approximation is not applicable at very low temperature (T <0.25 
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where Z(f3,ia) is obtained using eq. (|lj]) with a — > ia. The symmetric interval 
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MeV). So, we restrict ourselves to T > 0.25 MeV. Practically, for N = 10, the 
difference between p 9 pa /p and p 9 x /p is insignificant. On the other hand, for N=5 
the difference between p c Tpa /p and p c ex /p is noticeable at low T, whereas, this 
difference is small for iV=10. It implies that, though the grandcanonical partition 
function in SPA+RPA approach becomes almost exact, the small errors in its 
estimation causes a significant suppression of the level density obtained using 
number projected partition function (eq. at low T. Similar behaviour has 



been observed by Walt and Quick [|17| for the heat capacity. In other words, at 
low temperature, the small errors in grandcanonical partition function enhances 
the difference between Z N {(3) and Z N {(5). 

Let us now analyze the difference between canonical and grandcanonical treat- 
ments of level density and their deviations from the exact values. As discussed 
above (eqs. ^| and f|), in the present considerations, these treatments would differ 
only by the saddle point approximation to a single and a double Laplace-back trans- 
form. It is obvious from figures |l] and 0, the relative difference between the curves 
(a) and (b) or (c) and (d) at a fixed T and N is proportional to p c — p 9 . The dif- 
ference between the level densities using canonical and grandcanonical treatments 
decreases with increase in temperature or number of particles. For instance at T 
= 0.35 MeV, the ratio p c ex /p 9 x = 20.83 and 0.94 for N = 5 and 10, respectively, 
whereas, for T > 0.45 MeV this ratio becomes nearly equal to one. Finally, coming 
to the comparison between p and p, we find that p/p tends to unity with increase 
in number of particle and temperature. Thus, we can say that the saddle point 
approximation becomes almost exact with increase in temperature and particle 
numbers. 

In conclusion, we have analyzed the canonical and grandcanonical treatments 
of the nuclear level densities in SPA+RPA approach which accounts for large- 
amplitude thermal fluctuations and small-amplitude quantal fluctuations. The re- 



suits are compared with the corresponding exact ones. It is found that at low tem- 
perature the small errors in estimating the grandcanonical partition function cause 
significant suppression of level density obtained using number projected (canon- 
ical) partition function for small particle numbers. Also, we found that due to 
saddle point approximation to multiple Laplace-back transform, the grandcanon- 
ical treatment of the level density would be reliable at low temperature provided 
the particle number is relatively large. Furthermore, we must say that except at 
low temperatures where saddle point approximation itself is not reliable, a realis- 
tic value of level density can be obtained using SPA+RPA representation for the 
grandcanonical partition function. 
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Figure Captions 

Figure 1: Plot for p/p vs temperature for N = 5, where p represents the 
smoothed out exact (microcanonical) value of the level density and p corre- 
spond to the level density in a saddle point approximation which is obtained 
using (a) SPA+RPA approach to the canonical, (b) SPA+RPA approach to 
the grandcanonical, (c) exact canonical and (d) exact grandcanonical partition 
functions. 

Figure 2: Same as figure 1 but for N = 10. Note that the curves (a) and (b) 
or (c) and (d) are very close to each other at all temperatures. 
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